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SUMMARY 


A  large  experimental  study  has  been  undertaken  at  the  USAE  Waterways  Experiment 
Station  (WES)  as  part  of  the  ongoing  U.S.  Army  Engineer  Earthquake  Engineering 
Research  Program  (EQEN)  to  investigate  the  behavior  of  liquefying  soil  materials.  In 
this  study,  experiments  have  addressed  the  development  of  excess  pore  pressure  and 
the  onset  of  liquefaction  within  a  deep  soil  column  using  the  new  large  geotechnical 
centrifuge  and  earthquake  shaker.  The  approach  adopted  was  to  use  standard 
procedures  for  the  assessment  of  liquefaction  resistance  to  derive  values  of  the 
strength  reduction  factor  from  the  experimental  results,  with  the  object  of 
confirming  the  validity  of  the  simplified  Ka  approach  under  high  effective  overburden 
stresses.  Data  of  excess  pore  pressures  from  the  experiments  showed  the  soil  column 
reaching  a  state  of  initial  liquefaction  over  a  range  of  depths  up  to  around  20  m  (65  ft) 
under  moderate  levels  of  base  input  shaking,  and  from  this  data  values  of  could  be 
derived.  As  expected,  these  showed  considerable  variability,  being  sensitive  to  initial 
assumptions.  More  significantly,  at  greater  depths,  the  development  of  excess  pore 
pressure  was  capped  and  despite  continued  shaking  at  similar  amplitude,  the 
development  of  excess  pore  pressure  did  not  reach  a  sufficient  level  to  cause  initial 
liquefaction.  Similar  observations  have  been  found  in  the  literature  from  torsional 
shear  tests.  It  is  concluded  that  the  effects  of  confining  stress  and  the  strain  boundary 
conditions  which  exist  in  the  soil  column  in  the  field  are  significant  in  controlling  the 
development  of  excess  pore  pressure.  In  the  absence  of  initial  liquefaction,  the  use  of 
the  K<j  factor  to  assess  the  cyclic  resistance  of  the  soil  is  not  considered  appropriate. 
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1  INTRODUCTION 

The  current  state  of  practice  for  determining  the  performance  of  soils  that  undergo 
earthquake  induced  shear  strain  and  consequently  pore  water  pressure  build-up  in 
saturated  soils  is  to  determine  if  a  soil  will  or  will  not  liquefy.  The  performance  and 
safety  of  structures  are  based  on  rigid  sliding  block/slices  limit  equilibrium  methods 
to  determine  slip-plane  stability  and  deformation  using  residual  strengths  in  the  case 
of  liquefaction. 

In  the  standard  design  methods  for  the  assessment  of  liquefaction  resistance  of 
saturated  sands  under  high  effective  confining  stresses,  reliance  has  been  placed  on 
laboratory  test  data  which  show  a  general  reduction  in  the  cyclic  strength  of  the  soil. 
This  is  accounted  for  in  the  widely  used  'simplified  procedure'  through  the  use  of  the 
K<t  factor.  Seed  and  Harder  (1990).  At  high  effective  confining  stresses,  the  factor 
can  require  a  substantial  reduction  in  the  liquefaction  resistance  of  the  soil.  Studies  in 
the  literature  indicate  that  values  of  calculated  by  different  researchers  can  vary 
widely  and  that  extrapolation  of  the  simplified  procedure  to  depths  greater  than  15m  is 
not  supported  by  case  history  data  from  the  field. 

A  large  experimental  program  has  been  undertaken  at  the  USAE  Waterways 
Experiment  Station,  Vicksburg  Mississippi  using  the  new  centrifuge  research  facility 
and  earthquake  shaker. 
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2  THE  ARMY  CENTRIFUGE,  EARTHQUAKE  SHAKER  AND 
SPECIMEN  CONTAINER 

2. 1  The  Army  Centrifuge  and  earthquake  actuator 

The  design  specification  of  the  centrifuge  was  finalised  following  a  broad  ranging 
review  of  the  typical  field  structures  and  problems  with  which  the  Corps  is  principally 
engaged  Ledbetter  (1991).  The  design  of  the  centrifuge  was  based  on  the  French 
designed  Acutronic  661,  665,  and  680  series  of  geotechnical  centrifuges.  Similar  in  its 
physical  dimensions  to  the  680  model,  the  WES  centrifuge  (Fig.  2.1)  is  distinguished 
by  its  significantly  higher  performance  capability.  Key  characteristics  are 
summarised  in  Table  2.1. 


Figure  2. 1  The  Army  Centrifuge 


l™1 

Radius  to  platform 

6.5  m 

Payload  at  143  g 

8000  kg 

Payload  at  350  g 

2000  kg 

Capacity 

1144  g  tonnes 

Table  2.1  Key  performance  characteristics 

The  earthquake  actuator  used  in  the  present  series  of  experiments  was  a  large 
mechanical  shaker  designed  to  provide  a  single  frequency  input  motion  of  variable 
duration  to  the  base  of  the  specimen  container.  The  actuator  was  based  on  a  design 
successfully  developed  at  Cambridge  University  known  as  the  SAM  actuator, 
Madabushi  et  al.  (1996),  energy  is  stored  in  flywheels  in  the  form  of  angular 
momentum.  The  advantages  of  this  design  are  that  it  is  a  simple  mechanical  system 
capable  of  high  g  operation,  with  low  cost  of  construction.  Disadvantages  include  a 
restricted  ability  to  control  the  amplitude  of  input  motions  and  difficulties  in 
controlling  high  frequency  vibrations  under  low  amplitude  operations. 
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The  use  of  a  single  frequency  of  input  motion  was  considered  appropriate  to  the 
nature  of  the  experiments,  the  results  of  which  would  be  compared  with  laboratory 
cyclical  element  tests  and  standard  design  methods  for  liquefaction  assessment  which 
were  based  on  numbers  of  ‘equivalent’  uniform  cycles. 

The  general  arrangement  of  the  actuator  is  shown  in  Figs  2.2  and  2.3.  Fig  2.4  shows  a 
photograph  of  the  actuator  with  the  motor  on  the  right  and  high  vertical  support  walls 
enclosing  the  shaking  table  and  specimen  container  on  the  left. 

A  system  of  linkages  and  eccentrics  transfer  the  stored  energy  of  the  fly  wheels  to  the 
shaking  platform  and  thence  into  the  soil  specimen,  Figure  2.3.  A  hydraulic  or 
electrical  motor  drives  the  flywheels  up  to  full  speed,  and  then,  on  a  signal,  a  high 
speed  clutch  grabs  the  oscillating  shaft  and  transfers  energy  into  the  model  until 
another  signal  releases  it  again.  Clearly  the  frequency  of  the  oscillation  is  directly 
proportional  to  the  speed  of  the  motor  (and  flywheels).  The  amplitude  is  controlled  by 
the  arrangement  of  the  eccentrics;  three  displacement  amplitudes  for  the  platform  are 
available  (+/-  0.49mm,  +/-1.47mm  and  +/-4.41mm). 


-«<  ■ 


Figure  2.2  Elevation  of  the  earthquake  actuator,  Butler  (1999) 
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Figure  2.3  Plan  of  the  earthquake  actuator,  Butler  (1999) 


Figure  2.4  The  earthquake  actuator,  seen  from  the  rear 


The  shaker  is  designed  (structurally)  for  operation  up  to  150g,  at  which  the  maximum 
load  capacity  of  the  base  oscillating  platform  is  reached  (75  tonnes,  or  500  kg  at  lg). 
The  design  maximum  lateral  force  which  the  mechanism  may  exert  on  the  platform  is 
30  tonnes,  and  the  maximum  frequency  at  which  the  shaker  may  be  safely  operated  is 
150  Hz  (eg.  1Hz  prototype  at  150g,  or  3Hz  prototype  at  50g). 
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2.2  Specimen  container  (the  ESB  box) 

The  specimen  is  built  within  a  hollow  rectangular  model  container  (ESB)  comprising 
a  series  of  aluminium  alloy  rings  stacked  one  above  the  other,  and  separated  by  a 
elastic  medium,  Figure  2.5.  Several  of  these  chambers  have  been  constructed,  and 
extensive  dynamic  analysis  and  testing  has  been  carried  out  to  determine  their 
dynamic  response  characteristics,  Butler  (1999). 


Figure  2.5  The  WES  ESB  container,  Mk  II,  and  an  individual  aluminium  ring 


The  model  container  has  internal  dimensions  of  627mm  deep  by  315mm  wide  by 
796mm  long.  Each  of  the  eleven  aluminium  alloy  rings  is  50mm  high.  The  rings  are 
not  stiff  enough  along  their  long  dimension  to  support  the  outward  pressure  from  the 
soil  inside  under  high  g,  but  they  are  supported  by  the  massive  vertical  reaction  walls 
of  the  shaker  unit  itself.  A  rubber  sheet  separates  the  rings  from  the  steel  walls  on 
either  side.  This  concept  has  the  added  advantage  of  raising  the  centre  of  gravity  of 
the  reaction  mass  in  line  with  the  centre  of  gravity  of  the  specimen,  thus  minimising 
eccentric  forces  which  may  lead  to  rocking. 

Up  each  of  the  end  walls  of  the  ESB  are  positioned  thin  metal  sheets,  termed  shear 
sheets,  fixed  securely  to  the  base  of  the  chamber,  which  can  accommodate  the 
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complementary  shear  force  generated  by  the  horizontal  shaking  on  vertical  planes 
within  the  specimen  and  transmit  that  force  to  the  base  of  the  container,  Fig.  2.6.  This 
improves  the  uniformity  of  the  stress  field  at  each  elevation  along  the  model,  reducing 
the  tendency  for  the  chamber  to  ‘rock’. 


shear  sheet,  with 
rough  sand  finish, 
bolted  to  the  floor 
of  the  ESB  and 
laid  up  against  the 
end  wall 


stacked 

aluminium  rings, 
separated  by 
elastic  sealant 


Figure  2.6  Shear  sheets  at  each  end  of  the  ESB 

In  the  bottom  of  the  model  chamber,  slotted  metal  plates  were  used  to  create  flow 
channels  across  the  base  to  improve  the  uniformity  of  the  saturation  process.  These 
slots  were  filled  with  a  coarse  sand,  to  create  a  level  base  for  the  specimen.  To  ensure 
accuracy  of  weighing  the  specimen,  the  ports  and  slots  were  filled  with  fluid  prior  to 
the  dry  pluviation  of  the  specimen  proper.  The  combination  of  slots  and  coarse  sand 
was  designed  to  ensure  full  shear  transfer  into  the  specimen,  forming  a  ‘rough’  base, 
Fig.  2.7. 
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slotted  plate  in 
base  of  ESB  to 
assist  distribution 
of  pore  fluid  and 
provide  base 
roughness.  Under 
the  slots  can  be 
seen  the  metal 
porous  plate 


glass  sheets  form 
the  side  walls  of 
specimen 
container 


Figure  2.7  Slotted  plate  in  the  base  of  the  ESB 

The  ESB  concept  is  to  create  an  equivalent  shear  beam  with  an  average  stiffness 
comparable  to  the  stiffness  of  the  soil  specimen.  Expressed  rigorously,  the  concept  is 
more  accurately  defined  as  achieving  a  dynamic  response  which  does  not  significantly 
influence  the  behaviour  of  the  soil  specimen  inside.  For  certain  classes  of  experiment, 
it  may  be  expected  that  the  stiffness  of  the  soil  (at  least  near  the  end  walls)  would  not 
reduce  significantly  during  shaking  due  to  excess  pore  pressure  rise,  and  the  stiffness 
of  the  chamber  may  be  designed  accordingly,  perhaps  considering  a  shear  modulus 
appropriate  to  the  level  of  dynamic  strain  expected  in  the  soil  free-field  at  mid-depth. 
In  experiments  involving  the  liquefaction  of  large  volumes  of  soil  inside  the  container, 
the  stiffness  changes  (and  hence  dynamic  response)  will  change  dramatically 
throughout  the  base  shaking.  A  stiff  chamber  may  lead  to  undesirable  effects,  as 
observed  by  Peiris  (1999)  who  observed  that  liquefaction  in  a  loose  saturated  sand 
model  did  not  occur  near  the  stiff  end  walls  of  the  chamber.  A  chamber  with  no 
stiffness  simply  adds  mass  to  the  soil  specimen,  again  changing  its  dynamic  response. 
This  poses  a  particular  challenge  for  the  ESB  design. 

The  WES  ESB  was  assembled  using  a  urethane  adhesive  sealant  (commonly  used  as  a 
windshield  sealant  for  cars)  between  the  aluminium  alloy  rings,  a  material  which  was 
found  to  have  good  elastic  properties  and  which  bonded  well  to  the  metal  and  to  itself. 
The  new  container  has  a  low  shear  stiffness  and  a  first  mode  around  16Hz,  with 
second,  third  and  fourth  modes  at  48,  79  and  109Hz  respectively,  Butler  (1999).  [A 
typical  saturated  specimen  at  50g  in  the  WES  ESB  will  have  a  theoretical  natural 
frequency  around  92Hz,  based  on  the  small  strain  shear  modulus.]  The  stiffness  of 
the  ESB  was  then  investigated  to  assess  its  impact  on  the  resulting  combined  soil- 
container  response. 
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2. 3  Dynamic  response  of  the  soil  container  and  specimen 

In  his  doctoral  thesis,  Butler  (1999)  has  completed  a  thorough  theoretical  and 
experimental  analysis  of  the  dynamic  response  of  the  coupled  soil-container  system. 
At  high  g  the  soil  and  container  act  as  a  coupled  system,  where  the  lower  stiffness  of 
the  container  reduces  the  natural  frequency  (slightly)  of  the  combined  system 
compared  to  the  soil  column  alone.  However,  provided  the  driving  frequency  is  low 
relative  to  the  natural  frequency  of  the  coupled  system,  Butler  demonstrates  that  the 
displacement  response  of  the  system  is  unaffected  compared  to  the  soil  acting 
independently,  an  ideal  situation. 

For  higher  driving  frequencies,  Butler  concludes  that  it  would  be  necessary  to 
reconsider  the  elastic  stiffness  of  the  ESB  container,  and  to  tune  the  container  to 
ensure  that  even  with  the  expected  level  of  degradation  in  the  soil  specimen,  the 
coupled  system  did  not  deviate  significantly  from  the  condition  of  the  soil  column 
alone.  This  may  be  possible  by  adding  mass  to  the  rings  of  an  initially  stiffer  ESB  to 
reduce  its  first  mode  to  the  desired  level,  Butler  (1999). 

The  liquefaction  of  a  level  sand  bed  has  previously  been  the  subject  of  other  research. 
Experiments  were  conducted  at  many  centrifuge  centres  under  the  VELACS  project, 
Arulanandan  and  Scott  (1994).  The  objective  of  the  WES  study  is  to  investigate  the 
onset  of  liquefaction  under  much  higher  initial  effective  overburden  stresses. 
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Figure  2.8  Cross-section  through  Model  3c 

In  the  first  major  series  of  experiments  using  the  WES  earthquake  actuator,  the  ESB 
container  was  filled  with  saturated  sand,  forming  a  level  sand  bed,  Fig.  2.8.  These  are 
discussed  in  detail  below.  Instruments  were  placed  at  different  depths  to  measure  pore 
pressures  and  accelerations.  The  centrifuge  was  operated  at  50g,  creating  a  field 
equivalent  of  a  site  approximately  40m  long  by  15m  wide  by  over  26m  deep  (for  the 
deepest  specimens).  Surcharges  and  variations  in  the  depth  of  the  phreatic  surface 
have  been  used  to  achieve  higher  effective  overburden  stresses,  up  to  around  1000 
KPa. 

When  the  actuator  is  operated,  the  table  is  thrown  backwards  and  forwards  through  a 
fixed  amplitude  of  displacement,  which  is  governed  by  the  mechanical  gearing  of  the 
system,  as  described  above.  In  a  perfect  system,  the  peak  amplitude  would  be  reached 
in  the  first  half  cycle,  but  in  practice  tolerances  in  the  mechanical  design,  and  inertia 
of  the  specimen  and  container  interfere  with  the  theoretical  output  and  the  build-up  to 
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a  peak  motion  takes  several  cycles.  Higher  harmonics  of  the  driving  frequency  are 
also  often  present. 


Figure  2.9  Time  histories  from  Model  3c  Earthquake  1 

A  schematic  cross-section  through  a  typical  deep  (525mm)  model,  Model  3c,  is  given 
in  Figure  2.8.  Typical  earthquake  time  histories  from  Model  3c  are  illustrated  in 
Figure  2.9,  which  shows  the  base  shaking  motion,  motion  at  the  middle  of  the  loose 
layer,  motion  near  the  surface,  as  well  as  the  excess  pore  pressure  response  at  two 
depths  in  the  specimen  for  the  first  earthquake,  earthquake  1 . 

Butler’s  dynamic  analysis  of  the  container  and  specimen  in  his  doctoral  thesis  (1999), 
has  confirmed  that  the  soil-structure  system  is  responding  satisfactorily  without 
unacceptable  boundary  effects.  In  following  sections  of  this  paper,  the  data  of 
earthquake  response,  excess  pore  pressures  and  accelerations  are  analysed  in  the 
context  of  conventional  design  methods  for  liquefaction  assessment.  However,  the 
techniques  of  dynamic  analysis  applied  to  study  the  response  of  the  whole  system  may 
also  be  applied  with  benefit  to  individual  records  within  the  soil  specimen,  providing 
additional  valuable  insights  into  the  phenomenon  of  excess  pore  pressure 
development.  The  remainder  of  this  section  illustrates  how  this  technique  may  be  used 
to  study  the  onset  of  liquefaction  in  detail. 
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Figure  2.10  Harmonic  wavelet  plots  of  accelerometers  12610  (loose  layer)  Model  3c 
earthquake  1 
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Figure  2.1 1  Harmonic  wavelet  plots  of  accelerometers  12609  (base  of  sand)  Model 
3  c  earthquake  1 

Using  harmonic  wavelet  analysis,  time-frequency  ‘maps’  can  be  generated  from  the 
time  histories  of  acceleration  recorded  in  the  models.  Newland  and  Butler  have 
described  the  wavelet  approach  to  the  study  of  transient  signals  from  centrifuge 
models  in  recent  papers,  Newland  (1994,  1998),  Newland  and  Butler  (1998).  Figures 
2.10  and  2. 1 1  show  the  time-frequency  maps  for  ACC  126 10,  in  the  middle  of  the 
liquefying  loose  layer  compared  with  the  base  input  motion  ACC  12609  throughout 
the  earthquake. 
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Figure  2.12  Early  time  response  during  onset  of  liquefaction,  comparing  a  wavelet 
plot  and  time  histories  of  acceleration  with  pore  pressure.  Model  3c  earthquake  1 

During  the  early  stages  of  shaking,  there  is  a  rapid  increase  in  excess  pore  pressure 
and  consequent  reduction  in  soil  stiffness.  This  took  place  between  about  0.17  and 
0.27  seconds  in  earthquake  1 .  Time-frequency  maps  for  this  early  stage  is  shown  in 
Figure  2.12.  In  this  plot,  the  excess  pore  pressure  and  the  acceleration  in  the 
liquefying  layer  are  also  reproduced,  to  the  same  time  scale,  so  that  a  direct 
comparison  may  be  made.  It  is  clear  that  during  this  period  there  is  a  reduction  in  at 
least  one  of  the  main  resonance  frequencies  in  the  specimen,  from  around  90  to 
around  85  Hz.  This  frequency  corresponds  roughly  to  the  third  harmonic  of  the  input 


11 


Report  RSS/J9601R01  RevO 


frequency  (27  Hz)  and  also  to  the  likely  combined  soil  and  container  fundamental 
mode. 
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Figure  2.13  Differential  plot  of  comparison  between  accelerometer  12609  (positive) 
and  12610  (negative)  in  Model  3c  earthquake  1 

In  earthquake  1  and  in  later  earthquakes  a  sharp  rise  in  the  energy  at  around  140  Hz 
took  place  almost  exactly  coincident  with  excess  pore  pressure  reaching  100%.  This  is 
an  interesting  result  which  warrants  further  research.  140  Hz  is  near  the  fifth 
harmonic  of  the  input  shaking.  Examination  of  the  time  histories  in  Figure  2.9  showed 
sharp  spikes  developing  in  the  acceleration  record  as  the  soil  around  the  accelerometer 
liquefies.  The  physical  explanation  of  this  is  more  complex.  It  may  be  associated  with 
transient  shock  waves  ‘locking  up’  the  soil  instantaneously  and  then  ‘releasing’  it 
again.  If  the  shock  wave  was  generated  by  the  fundamental  shaking  motion  of  the 
specimen,  then  this  would  explain  why  the  high  frequency  spike  only  occurs  on  each 
half  cycle  of  base  motion.  Alternatively,  it  may  be  associated  with  a  fracture  type 
phenomenon  as  the  soil  liquefies  and  high  frequency  energy  is  generated. 

Additional  insights  may  be  gained  by  subtracting  the  absolute  values  of  two  harmonic 
wavelet  transforms  to  give  a  differential  plot.  This  gives  a  comparison  of  the  energy  at 
different  locations  in  the  model.  Figure  2.13  shows  a  differential  plot  for 
accelerometers  12609  and  12610  in  Model  3c,  earthquake  1.  The  intense  band  along 
the  base  of  the  figure  shows  that  the  input  shaking  energy  is  much  greater  at  the  base 
of  the  model  than  in  the  middle  of  the  loose  layer  (ACC  12610),  as  expected. 
However,  energy  at  around  84Hz  was  amplified  between  the  base  of  the  sand  bed  and 
the  middle  of  the  loose  layer,  a  small  distance  above.  Here  the  absolute  value  of  the 
differential  is  negative.  This  effect  was  transient  in  earthquake  1 ,  decaying  with  time, 
but  this  pattern  was  different  in  subsequent  earthquakes. 
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3  MODEL  TEST  CONFIGURATION  AND  INSTRUMENTATION 

3. 1  Outline  of  experiments 

The  model  tests  were  designed  to  investigate  the  liquefaction  of  a  loose  saturated 
layer  under  varying  effective  overburden  pressures.  The  principle  aim  of  the 
experiments  was  to  achieve  an  improved  understanding  of  the  factor  in 
liquefaction  analysis  through  centrifuge  model  tests  of  a  level,  saturated  sand  bed 
under  strong  base  shaking.  The  objective  of  the  series  was  to  capture  data  of 
accelerations  and  excess  pore  pressures  in  a  loose  layer  as  excess  pore  pressures  reach 
a  condition  of  initial  liquefaction  under  a  range  of  different  initial  effective 
overburden  stresses  ranging  from  1  tsf  to  10  tsf.  The  experiments  were  conducted  in 
the  Equivalent  Shear  Beam  (ESB)  model  container  using  the  new  earthquake  actuator 
on  the  WES  centrifuge,  described  above. 

3. 2  Summary  of  Model  Test  Series 

Table  3.1  summarises  the  experiments  conducted  during  1998  and  1999.  The  models 
are  grouped  in  series,  where  each  series  corresponds  to  a  different  target  range  of 
vertical  effective  overburden  stress  in  the  loose  layer.  In  all  cases,  the  bottom  1 60mm 
of  the  specimen  was  ‘loose’  and  the  upper  portion  ‘medium-dense’.  All  models  were 
shaken  at  50g.  Some  models  were  overconsolidated  by  a  factor  of  2.5  prior  to  shaking 
(achieved  by  running  the  centrifuge  upto  125g). 


Model 

series 

Models  in  series 

Effective 
overburden 
stress  in  loose 
layer 

Depth  of 
prototype 
(approx) 

Depth  of 
specimen 

Notes 

1 

a 

1  tsf 

15  m 

300  mm 

Ottawa  sand 

2 

a,  b,  c,  d,  e,  f 

1  tsf 

15  m 

300  mm 

Nevada  sand 

3 

a,  b,  c,  d,  e 

2  tsf 

26  m 

525  mm 

Nevada  sand 

4 

a,  b,  c,  d 

3  -  5  tsf 

26-40  m 

525  mm 

Nevada  sand 
with  lowered 
w.t.  or 
surcharge 

5 

a,  b,  c,  d 

7-  10  tsf 

54  -  63  m 

525  mm 

Nevada  sand 
with  lead 
surcharge 

Table  3.1  Summary  of  Model  series 

Table  3.2  provides  further  details  on  each  experiment,  including  the  densities 
achieved  in  the  specimen,  the  date,  the  number  of  earthquake  events  and  the  actual 
initial  vertical  effective  stress  in  the  centre  of  the  loose  layer.  Each  centrifuge 
experiment  generated  a  large  body  of  data  from  instrumental  records,  from 
measurements  made  during  construction  and  excavation,  from  photographs,  video 
records  and  other  sources. 
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Model 

Code 

WES 

CRC 

Ref. 

Overall 

Depth 

(mm) 

Density 

ov’ 

(tsf) 

Date 

OCR 

Earth¬ 

quakes 

Comments 

la 

1 

300 

52%  loose 
70%  dense 

1 

4/12/97 

1 

2 

Ottawa  sand,  Cambridge 
ESB 

2a 

2 

300 

44%  loose 
83%  dense 

1 

20/2/98 

1 

3 

Nevada  sand,  Cambridge 
ESB 

2b 

5 

300 

50%  loose 
75%  dense 

1 

30/4/98 

1 

2 

Nevada  sand,  ESB  #1 

2c 

8 

300 

49%  loose 
74%  dense 

1 

6/6/98 

1 

5 

Nevada  sand,  ESB  #1 

2d 

9 

300 

50%  loose 
75%  dense 

1 

12/8/98 

1 

4 

Nevada  sand,  ESB  #2 

2e 

10 

300 

49%  loose 
73%  dense 

1 

26/8/98 

2.5 

4 

Nevada  sand,  ESB  #1 

2f 

14 

300 

50%  loose 
75%  dense 

1 

22/9/98 

2.5 

4 

Nevada  sand,  ESB  #1 

3a 

3 

525 

34%  loose 
73%  dense 

2 

30/3/98 

1 

2 

Nevada  sand,  Cambridge 
ESB 

3b 

6 

525 

49%  loose 
77%  dense 

2 

1/5/98 

1 

3 

Nevada  sand,  ESB  #2,  no 
data  recovered  from  eql 

3c 

7 

525 

49%  loose 
79%  dense 

2 

5/5/98 

1 

3 

Nevada  sand,  ESB  #2 

3d 

3e 

11 

21 

525 

54%  loose 
80%  dense 
50%  loose 
75%  dense 

2 

3/9/98 

2.5 

1 

4 

Nevada  sand,  ESB  #2 

No  useful  data  recovered 

4a 

12 

525 

49%  loose 
80%  dense 

3 

14/9/98 

1 

4 

Saturated  to  top  of  loose 
layer  only.  Nevada  sand, 
ESB  #1 

4b 

13 

525 

56%  loose 
74%  dense 

3 

15/9/98 

2.5 

4 

Saturated  to  top  of  loose 
layer  only.  Nevada  sand, 
ESB  #2 

4c 

17 

525 

50%  loose 
75%  dense 

4.7 

22/10/98 

1 

4 

Lead  surcharge,  Nevada 
sand,  ESB  #1 

4d 

18 

525 

50%  loose 
68%  dense 

4.7 

4/11/98 

2.5 

4 

Lead  surcharge,  Nevada 
sand,  ESB  #1,  model 
subject  to  further 
earthquakes  to  test  new 
motor  prior  to  dissembly 

5a 

15 

525 

5 1  %  loose 
72%  dense 

7.4 

7/10/98 

1 

4 

Lead  surcharge,  Nevada 
sand,  ESB  #1 

5b 

16 

525 

49%  loose 
76%  dense 

7.4 

14/10/98 

2.5 

4 

Lead  surcharge,  Nevada 
sand,  ESB  #1 

5c 

19 

525 

52%  loose 
75%  dense 

9.2 

11/12/98 

1 

3 

Lead  surcharge,  Nevada 
sand,  ESB  #2 

5d 

20 

525 

57%  loose 
80%  dense 

9.2 

17/2/99 

1 

1 

Lead  surcharge,  Nevada 
sand,  ESB  #2 

Table  3.2  Summary  of  experiments 


The  objective  for  the  upper  portion  of  the  specimen  was  a  create  a  medium  dense  (Dr 
=  75%)  sand,  and  for  the  lower  layer  a  target  relative  density  of  around  50%.  The 
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depth  of  the  upper  portion,  water  table  and  magnitude  of  surcharge  were  calculated  to 
ensure  the  target  effective  overburden  pressure  was  achieved  at  mid-height  in  the 
loose  layer.  The  first  model  was  constructed  with  Ottawa  sand  (Model  la)  but  all 
subsequent  models  (which  included  repeats  of  earlier  experiments)  used  Nevada  sand 
(Model  series  2,  3  etc). 


796 


796 
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160 
r  15 


water  table  at 
soil  surface 


medium  dense  sand  layer, 
,  RD  =  75%. . : 


365 


160 

15 


Typical  cross  section  through 
Model  series  1  and  2 
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;;  V  ?> 
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-  t  r 
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medium  dense  sand  layer, 
,  RD  =  75%  . 


365 


160 
:  15 


Typical  cross  section  through  Models  4a,  4b  Typical  cross  section  through  Models  4c,  4d  and  5 


The  importance  of  the  target  elevation  in  the  middle  of  the  loose  layer  was  that  a  fixed 
volume  of  sand  of  similar  relative  density  lay  below. 

The  specimens  were  all  300mm  wide,  between  the  vertical  side  walls  of  the  ESB.  The 
side  walls  comprised  smooth  glass  plates  mounted  against  the  aluminium  rings, 
providing  a  sheer  low  friction  interface.  (This  was  in  contrast  to  the  end  walls,  which 
comprised  roughened  steel  ‘shear  sheets’,  described  above.)  On  the  base,  slotted 
plates  were  used  to  assist  the  saturation  process  by  distributing  fluid  over  the  base  of 
the  model.  Coarse  sand  was  used  to  fill  the  slots,  providing  a  filter  and  creating  high 
friction  interface  between  the  specimen  and  the  base  of  the  container. 

As  the  prototype  depth  was  increased,  the  options  were  to  shake  the  models  at  a 
higher  g  level,  to  depress  the  water  table  or  phreatic  surface,  or  to  use  surcharging  to 
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achieve  the  target  overburden.  It  was  determined  that  the  models  should  all  be  shaken 
at  50g,  and  that  surcharging  would  be  adopted  to  reach  the  required  stress  levels.  In 
the  Model  4  series,  two  models  used  a  depressed  phreatic  surface  and  two  used 
surcharging.  In  the  Model  5  series,  the  prototype  depth  was  too  great  to  be  achieved 
without  surcharging.  The  surcharge  comprised  lead  strips  laid  lengthwise  along  the 
surface  of  the  specimen,  with  additional  lead  plate  on  top  where  necessary. 

3.3  Materials 

The  materials  used  in  the  model  were  characterised  by  standard  laboratory  tests  to 
determine  parameters  such  as  dry  density  and  gradation.  With  the  exception  of  Model 
la,  all  models  were  constructed  using  Nevada  sand  and  saturated  with  a  glycerine- 
water  solution.  Tables  3.3,  3.4  and  3.5  present  key  material  parameters  for  Ottawa 
sand,  Nevada  sand  and  glyerine-water  solution  respectively. 


Specific  gravity 

2.68 

Maximum  void  ratio 

0.7633 

Minimum  void  ratio 

0.4762 

D50 

0.12  mm  (approx) 

Dio 

0.075  mm  (approx) 

Table  3.3  Ottawa  Sand  specification  (from  VELACS) 


Specific  gravity 

2.64 

Maximum  void  ratio 

0.757  (density  93.8  pcf) 

Minimum  void  ratio 

0.516  (density  108.7  pcf) 

D50 

0.18  mm  (approx) 

Dio 

0.1 1  mm  (approx) 

Table  3.4  Nevada  Sand  specification  (as  measured) 


The  pore  fluid  comprised  a  mixture  of  glycerine  and  water,  80%  by  weightfor 
experiments  conducted  at  50g.  Measurements  of  the  viscosity  of  glycerine-water 
mixes  at  a  range  of  temperatures  and  proportions  show  that  the  viscosity  is  sensitive 
to  both  parameters,  Fig  3.1,  Steedman  (1999). 
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Figure  3.1 

The  density  of  a  glycerine-water  mix  was  calculated  from: 


pm  =  pg(mg  +  mw)/(  m„  +  pg  mw) 


where  pm  is  the  density  of  the  mix,  pg  is  the  density  of  glycerine,  mg  is  the  mass  of 
glycerine,  and  mw  is  the  mass  of  water.  Table  3.5  summarises  the  properties  of  the 
glycerine-water  solution  used  as  the  pore  fluid. 


Density 

1200  kg/m3 

Viscosity 

50  cs 

Specific  Gravity 

1.26 

Composition 

80%  glycerine-water  mix  (by  weight) 

Table  3.5  Specification  for  pore  fluid  (as  measured) 


The  models  were  poured  dry  from  a  hopper  and  saturated  under  vacuum,  or 
occasionally  under  gravity.  Instrumentation  was  placed  in  the  model  as  it  was  being 
constructed. 

3. 4  Instrumentation  layout 

Each  specimen  was  instrumented  with  accelerometers  and  pore  pressure  transducers. 
Details  of  the  exact  locations  of  all  instruments  are  contained  in  the  factual  report 
prepared  for  each  experiment.  Figs.  3.2, 3.3,  and  3.4  show  the  typical  locations  and 
type  of  instrumentation  in  three  models,  representative  of  the  Model  2,  Model  3  and 
Model  5  series. 
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Figure  3.2  Instrumentation  typical  of  Model  2  series 
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Figure  3.3  Instrumentation  layout  typical  of  Model  3  series 
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Figure  3.4  Instrumentation  typical  of  Model  series  4  and  5 

The  accelerometers  used  in  the  models  were  piezoelectric  capacitative  devices 
manufactured  by  D  J  Birchall  (UK)  and  Endevco  (US).  In  the  above  figures,  the 
Birchall  devices  are  identified  by  a  four  digit  reference  number;  the  Endevco  devices 
are  numbered  12xxx.  Typical  calibrations  for  the  Birchall  devices  are  around  0.12 
mV/g.  The  Endevco  accelerometers  are  around  half  of  this  sensitivity,  with 
calibrations  typically  around  0.05  mV/g.  Cables  from  the  devices  were  routed  to  the 
side  walls  of  the  model  chamber  and  up  the  glass  sides,  Fig.  3.5. 
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Figure  3.5  Polaroid  photo  of  instrumentation  placement  at  layer  4  in  Model  4d 

Fig.  3.5  also  shows  clearly  the  side  wall  arrangement  used  in  the  experiments,  with 
glass  sheets  mounted  inside  the  aluminium  rings  of  the  ESB.  A  soft  elastic  sealing 
compound  was  used  to  fill  voids  between  the  glass  and  the  aluminium,  shown  in  black 
in  the  photograph. 

Vibrations  of  the  cables  was  minimised  by  bundling  them  together  and  securing  them 
to  gantries.  Amplifiers  mounted  in  a  junction  box  on  the  vertical  support  wall  of  the 
shaker  converted  the  input  signals  to  voltages. 

The  pore  pressure  transducers  used  in  the  experiments  were  manufactured  by  Druck, 
and  are  widely  used  in  centrifuge  modelling.  Typical  calibrations  depend  on  the  range 
of  the  device,  and  vary  from  around  0.17  mV/V/psi  for  a  5  bar  device  to  0.036 
mV/V /psi  for  a  35  bar  ppt.  Using  amplifiers  in  the  junction  boxes,  these  signals  were 
conditioned  and  stored  locally  on  an  on-board  computer  before  transmitting  them 
through  the  slip  rings  to  the  data  logging  system  in  the  control  room. 
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4  ANALYSIS  AND  INTERPRETATION  OF  EARTHQUAKE  MODEL 
TEST  DATA 

4. 1  State-of-practice  for  liquefaction  assessment 

The  established  method  of  assessment  of  liquefaction  resistance  of  a  soil  layer  or 
column  is  to  conduct  field  (in-situ)  investigations  and  to  correlate  measurements  of 
standard  penetration  resistance,  cone  resistance  or  other  parameters  with  field 
observations.  Estimation  of  the  cyclic  stress  ratio  (CSR)  and  the  cyclic  resistance 
ratio  (CRR)  are  then  used  to  assess  the  likelihood  of  liquefaction.  These  are  defined  as 
the  seismic  demand  on  the  soil  and  the  capacity  of  the  soil  to  resist  liquefaction 
respectively,  Youd  and  Idriss  (1997).  Assessment  of  the  CRR  is  subject  to  corrections 
such  as  those  proposed  by  Seed  (1983)  to  account  for  factors  such  as  depth  (Kct)  and 
high  initial  static  shear  stress  (Kct)  where  the  method  is  extended  beyond  the  original 
data-set.  Because  of  the  absence  of  case  history  data  against  which  comparisons 
could  be  made,  these  factors  have  had  to  be  developed  using  laboratory  test  data. 

The  present  state-of-practice  for  determining  the  performance  of  soils  that  undergo 
earthquake-induced  shear-strain  and  consequently  pore  water  pressure  buildup  in 
saturated  soils  is  to  determine  if  a  soil  will  or  will  not  liquefy.  The  performance  and 
safety  of  structures  are  based  on  rigid  sliding-block/slices  limit-equilibrium  methods 
to  determine  slip-plane  stability  and  deformation  using  residual  strengths  in  the  case 
of  liquefaction. 

In-situ  methods  such  as  Standard  Penetration  Test  (SPT),  Cone  Penetration  Test  and 
shear  wave  velocity  measurements  are  used  to  define  the  triggering  of  liquefaction. 
The  potential  for  liquefaction  relies  on  empirical  correlation  between  the  penetration 
resistance  and  the  performance  of  soil  deposits  in  past  earthquakes  (Seed,  1979).  The 
database  on  triggering  liquefaction  is  based  on  data  from  level  ground  where  surface 
evidence  occurred.  Soil  conditions  were  shallow  at  overburden  pressures  less  than 
about  95.8  kPa  (about  4.6m  of  depth).  The  data  base  is  normalized  to  an  earthquake 
of  magnitude  7.5.  Corrections  for  other  earthquake  magnitudes  have  to  be  made 
(Seed  &  Harder,  1990,  NCEER,  1996). 

To  extend  the  database  to  depths  representative  of  foundation  soils  under  dams 
requires  the  use  of  correction  factors  (Seed  &  Harder,  1990)  for  the  effects  of 
overburden  pressure  and  initial  static  shear-stress.  Correction  for  high  overburden, 

Ka,  is  based  on  laboratory  test  results  of  the  ratio  of  cyclic-shear-stress,  CSS,  to  cause 
liquefaction  at  an  overburden  effective-stress  state,  a'0,  to  that  at  a  o'0  =  95.8  kPa. 
Correction  for  initial  static  shear-stress,  Ka,  is  based  on  the  ratio  of  CSS  to  cause 
liquefaction  with  initial  static  shear-stress  applied  to  that  at  no  initial  static  shear- 
stress. 

A  factor  of  safety,  FSL,  against  the  occurrence  of  liquefaction,  defined  as  100  percent 
pore  water  pressure,  can  be  calculated  as: 

FSL  =  f  [(N,)60]  Km  Ka  Ka  /xeq 

Where  (Ni)6o  is  from  SPT  relations.  Km  is  earthquake  magnitude  or  duration 
correction  factor  and  xeq  is  a  percentage  of  peak  earthquake  induced  shear  stress. 
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The  FSL  can  be  related  to  percent  excess  pore  water  pressure. 

The  correction  factor  Ka  has  a  large  influence  for  dam  foundations.  It  can  reduce  the 
CSS  to  cause  liquefaction  to  about  45  percent  of  its  in-situ  value  at  pressures  about 
670  KPa  that  would  exist  beneath  an  embankment  dam  about  30.5m  high. 
Alternatively,  if  SPT  criteria  were  being  derived  for  remedial  treatment  of  this  dam  to 
limit  the  potential  level  of  earthquake  generated  pore  water  pressures,  the  Ka  factor 
would  cause  an  increase  of  more  than  double  in  the  required  penetration  values  to  be 
measured  in  the  field.  Clearly,  the  correction  factor  can  have  a  major  impact  on  the 
potential  for  triggering  liquefaction  or  the  excess  residual  pore  water  pressures  and  on 
the  cost  of  remediation. 

The  correction  factor  Ka  can  also  contribute  significantly  to  the  reduction  of  the  in- 
situ  strength.  However,  for  relative  densities  above  45-50  percent,  Ka  can  have  a 
positive  effect  on  the  in-situ  strength.  The  Ka  relationships  are  not  well  defined. 

In  application  of  the  state-of-practice,  the  following  assumptions  are  made  inherently: 
(1)  the  soil  is  always  undrained,  (2)  liquefaction  occurs  instantaneously  and  the  soil 
shear  strength  jumps  to  residual  state,  (3)  residual  strength  is  constant  with  monotonic 
loading,  (4)  liquefaction  is  independent  of  soil  zone  thickness,  permeability,  or 
boundary  conditions,  (5)  liquefaction  is  independent  of  when  the  earthquake  peak 
energy  arrives,  (6)  behavior  of  the  liquefied  soil  and  its  resultant  effects  on  a  dam  are 
independent  of  the  soil  zone  thickness,  permeability,  and  boundary  conditions,  (7) 
dam  stability  and  deformation  are  controlled  by  slip-planes  independent  of  the 
liquefied  soil  zone  thickness  and  behavior,  and  (8)  non-liquefied  soil  at  a  site  is 
unaffected  by  the  earthquake.  Field  behavior,  numerical  analyses,  and  physical  model 
tests  show  that  these  assumptions  are  invalid. 

4.2  Research  problem 

Current  studies  for  more  thorough  evaluation  of  liquefaction  and  for  remediation 
design  and  analysis  have  shown  serious  limitations  in  the  state-of-practice.  The  state- 
of-practice  can  force  costly  excessive  remediations  to  be  undertaken  when  possibly  no 
action  is  required,  but  it  can  also  lead  to  unsafe  conclusions  in  other  cases. 

Significant  progress  is  being  made  in  the  development  of  numerical  methods  for 
analysis  of  liquefaction  and  the  consequences.  However,  the  engineering  profession 
will  most  likely  always  use  empirical  correlation  (Seed’s  or  others)  of  in-situ 
measurements  versus  potential  liquefaction,  pore  water  pressure  generation  and 
earthquake  response.  Every  time  a  site  is  evaluated  for  a  seismic  design  or  potentially 
liquefiable  soil  is  improved  and  a  dam  remediated,  in-situ  measurements  will  be 
made.  A  value/range  of  in-situ  measurement  to  achieve  will  be  specified  for  a 
construction/remediation  contractor.  Some  in-situ  measure  will  be  used  to  judge  soil 
conditions/improvement  and  seismic  safety  of  a  dam  or  site.  Therefore,  improvement 
in  the  current  state-of-practice  and  the  empirical  correlations  between  in-situ  measures 
and  performance  of  soil  deposits  has  to  be  made. 

Current  studies  of  the  seismically-induced  deformation  behavior  of  dams  indicate  that 
as  soils  are  progressing  toward  liquefaction  (pore  water  pressure  is  increasing  and 
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shear  strain  is  occurring)  significant  deformations  of  the  structure  can  occur.  Failure 
(damaging  levels  of  deformation)  can  develop  significantly  before  the  complete  initial 
liquefaction  stage  (100  percent  pore  water  pressure  ratio)  is  reached.  Depending  on 
specific  conditions  involving  the  location,  depth  and  extent  of  liquefying  soil  and  the 
driving  forces,  a  structure  may  fail  at  only  50  percent  strength  reduction.  In  this  case, 
remediation  to  assure  safe  performance  is  required  to  prevent  serious  damage 
significantly  before  an  initial  liquefaction  condition  and  a  residual  strength  stage  are 
arrived  at  in  the  soil. 

The  problems  in  the  current  state-of-practice  stem  mainly  from  the  fact  that  for  the 
sites  that  have  liquefied  and  constitute  the  empirical  basis  for  analysis,  the  following 
are  not  known:  (1)  the  exact  and  complete  soil  conditions  and  profiles,  (2)  the  real 
behavior  that  occurred  in  the  soils  during  and  after  the  earthquakes  or  the  various 
influences  on  the  behavior,  (3)  the  development  of  pore  pressures  or  strains  in  the 
assumed  non-liquefied  soils  (used  in  comparison  to  liquefied  soils)  which  may  have 
changed  state  during  the  earthquakes,  and  (4)  whether  artificial  and  possibly  incorrect 
conditions  in  laboratory  testing  may  have  led  to  conclusions  not  totally  applicable  to 
the  field  behavior. 

Improved  definition  and  physical  evidence  is  therefore  needed  of  the  processes  and 
mechanisms  involved  as  a  soil  progresses  to  liquefaction  and  residual  strength.  This 
is  needed  to  allow  refined  analyses  for  dam  safety  and  more  cost-effective  and  safe 
remediation  design  and  analysis.  Because  various  assumptions  can  be  made  coupled 
with  methodologies  and  numerical  analyses  that  can  give  solutions  or  answers  to 
almost  anything,  the  reality-check  of  solutions  must  come  from  field  or  equivalent- 
field  data  of  behavior  under  well  known  and  defined  conditions. 

The  earthquake  response  database  needs  to  be  expanded  with  more  complete  data  to 
provide:  (1)  the  necessary  advance  in  the  state-of-practice,  (2)  a  basis  for  modification 
and  improvement  of  current  methodology  and  assumptions,  and  (3)  definition  of  the 
physical  processes  and  mechanisms  involved  in  the  liquefaction  process  and  resultant 
effects  on  dam  behavior.  This  would  also  provide  the  fundamentals  and  basis  for 
development  of  new  methodology  and  analyses.  New  methodologies  have  to  be 
based  on  correct  mechanisms  and  processes. 

Current  specific  needs  for  more  thorough  earthquake  engineering  analyses  can  be 
identified  from  examination  of  the  last  two  decades  of  experience  in  seismic 
evaluation  of  embankment  dams  and  the  serious  limitations  that  arise  when 
remediation  design  is  attempted.  Studies  involving  liquefaction,  stability,  and 
seismic-induced  deformation  behavior  of  dams  raise  serious  questions  that  impact  the 
safe  performance  and  needs  for  remediation  (e.g.,  Ledbetter  and  Finn,  1993, 

Ledbetter,  et.  al.,  1994,  Finn  and  Ledbetter,  1991,  Finn,  et.  al.,  1991  and  1994,  Vaid 
and  Chem,  1985,  Vaid  and  Thomas,  1994,  Bryne  and  Harder,  1991).  For  example, 
Vaid  and  Thomas  (1994)  demonstrated  that  Ko  for  specific  sand  types  may  be 
substantially  less  (more  than  a  factor  of  2)  than  suggested  by  Seed  and  Harder.  Pillai 
and  Byrne  (1994)  showed  for  the  foundation  materials  of  Duncan  Dam,  a  Kcr  of  0.6 
compared  with  0.4  from  the  Seed  and  Harder  relation.  This  made  a  difference  between 
recommending  remediation  and  no  remediation.  There  is  a  large  spread  of  Ka 
relations  and  data  by  various  researchers  that  can  influence  whether  to  remediate  a 
dam  at  costs  in  the  tens  of  millions  of  dollars  or  not  remediate. 
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For  the  past  thirty  years,  research  has  primarily  concentrated  on  the  triggering  of 
liquefaction  of  soils  both  in  the  laboratory  and  in  field  sites.  A  significant  database 
has  resulted  of  very  important  and  necessary  information  concerning  the  stress-state 
triggering  of  liquefaction,  cyclic  load  stress-ratio  and  the  dynamic  properties  of  soils. 
Newmark,  1965,  stressed  that  what  counted  was  whether  the  deformations  that  a  dam 
suffered  during  an  earthquake  were  tolerable  or  not.  Peck,  1992,  stated  that  of  all 
measures  of  safety,  the  most  directly  applicable  results  are  the  anticipated 
deformations. 

Due  to  the  lack  of  knowledge  and  experience  involving  the  behavior  of  liquefiable 
soils  under  field  conditions,  simplifying  assumptions  are  forced  to  be  made  for  the 
critical  safety  of  dams,  concerning  behavior.  Possible  significant  controlling 
influences  such  as  permeability,  boundary  layers  and  pore  water  pressure  dissipation 
and  migration  must  be  disregarded. 

Some  specific  needs  are:  (1)  well  defined  and  complete  shear  stress-strain  response 
curves  for  earthquake  loading  including  the  residual  strength  portion,  (2)  strains 
within  a  problem  soil  mass,  (3)  effects  of  soil  zone  thickness,  permeability,  and 
boundary  conditions,  (4)  influence  of  adjacent  soil  materials  and  of  their 
permeabilities,  (5)  dissipation  and  movement  of  excess  residual  pore  pressure  both 
during  and  after  an  earthquake,  (6)  redistribution  of  stresses  as  a  soil  is  losing 
strength,  (7)  interaction  of  remediation  materials  and  adjacent  soil,  (8)  dynamic 
response  of  remediation  materials  and  of  remediated  zones,  (9)  improved  Ka  and  Ka 
factors  for  the  field  evaluation  of  remediation  achievement  and  for  improved  first 
estimates  of  liquefaction  potential,  (10)  effects  of  strong  aftershocks  and  (11)  dam 
internal  behavior  and  failure  mechanisms  in  response  to  earthquake  loading  and 
strength  degradation. 

The  experimental  investigation  using  the  Army  centrifuge  at  WES  has  been  used  to 
address  these  issues  as  it  has  provided,  for  the  first  time,  a  substantial  body  of  data  of 
the  behavior  of  liquefying  soil  deposits  under  a  wide  range  of  initial  effective 
confining  stresses.  In  the  following  sections,  the  method  by  which  the  data  has  been 
analysed  is  presented.  Although  the  experiments  adopted  a  target  effective  confining 
stress  within  a  loose  layer  at  some  depth,  instrumentation  at  other  depths  also 
provided  valuable  information  on  the  degradation  of  the  soil  column.  In  some  cases, 
upper  layers  were  clearly  affected  by  the  generation  of  excess  pore  pressure  in  layers 
below.  This  is  typical  of  field  behavior,  although  not  of  laboratory  element  tests.  The 
first  stage  in  the  experiment  analysis  was  to  develop  the  method  by  which  Ka  could 
be  calculated  within  the  soil  column. 

4.3  Calculation  of  Ka  from  measurements  made  in  a  soil  column 

The  stress  history  at  the  location  of  each  transducer  was  calculated  using  a  time 
history  of  acceleration.  One  accelerometer  record  was  selected  based  on  the  quality  of 
its  signal  and  its  location.  Generally  the  device  was  located  on  the  bottom  ring  of  the 
ESB,  but  occasionally  the  accelerometer  at  the  base  of  the  sand  specimen  was  used 
instead.  In  some  cases,  where  amplification  was  considered  significant,  an 
amplification  analysis  was  carried  out  by  examining  the  peak  acceleration  at  different 
elevations  in  the  model,  based  on  the  very  early  cycles  of  shaking  before  significant 
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degradation  had  taken  place.  The  mean  of  positive  and  negative  cycles  was  used,  and 
commonly  high  frequency  spikes  were  removed  by  smoothing  where  these  were 
considered  to  be  unrepresentative  of  the  character  of  the  motion.  Plotting 
amplification  as  a  function  of  the  base  input  motion  revealed  that  generally  the 
amplification  was  around  1,  Fig  4.1. 
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Note:  amplification  calculated  at  each  elevation  as  the  mean  of  the  positive 
and  negative  peak  values  during  the  second  cycle  of  shaking,  betw  een 
0.28  and  0.32  seconds.  The  mean  of  both  accelerometers  at  each 
elevation  was  then  used,  with  the  exception  of  ACC  12610  (unreliable) 
and  on  the  base,  w  here  only  one  accelerometer  w  as  available. 


Note:  amplification  calculated  at  each  elevation  as  the  mean  of  the  positive 
and  negative  peak  values  during  the  second  cycle  of  shaking,  between 
0.95  and  1.1  seconds.  The  mean  of  both  accelerometers  at  each  elevation 
was  then  used,  w  ith  the  exception  of  the  near  surface,  the  middle  of  the 
loose  layer  and  on  the  base,  w  here  only  one  accelerometer  w  as  avaSable. 


Model  2b,  amplification  during  second  cycle  of  base  shaking 


Model  2c,  amplification  during  early  cycles  of  shaking 


Figure  4.1  Amplification  during  early  cycles  of  shaking  for  two  models 

The  cyclic  stress  ratio  (CSR)  at  any  depth  was  calculated  using  the  relationship 
originally  proposed  by  Seed  and  Idriss  (1971)  and  developed  by  Ishihara  (1996), 


CSR  =  =  0.65  rH  ^ 

f  a  » 

°v  g 

where  amax  is  the  peak  acceleration  at  the  ground  surface  and  xav  is  the  shear  stress  on 
a  horizontal  plane. 
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Comparison  between  'effective' 
stress  ratios  calculated  for  deep 
soil  columns  using  data  from 
SHAKE-7  analyses  of  two  field 
earthquakes  (Loma  Prieta  and  El 
Centro),  one  artificial  earthquake 
(Folsom  Record  B)  and  one 
centrifuge  earthquake  (Model  3c 
earthquake  1).  Soil  column  was 
assumed  to  be  all  sand,  shear 
modulus  and  damping  properties 
set  to  'average'  values.  The  water 
table  assumed  to  be  at  the  surface. 
Graphs  also  show  'effective'  stress 
ratio  calculated  using  ground 
surface  acceleration  based  on  two 
alternative  reduction  factors  with 
depth, 

rd  =  l-0.015z,  and 
rd  =1.25 -0.2  In  (z) 


Figure  4.2  Output  from  SHAKE  computations  of  the  response  of  a  soil  column  to 
four  different  earthquake  motions 

There  is  a  wide  range  of  reduction  factors  rd  proposed  in  the  literature.  For  this  study, 
a  mean  rd  was  developed  based  on  analyses  using  SHAKE  of  soil  columns  similar  to 
those  constructed  in  the  centrifuge  experiments.  Figure  4.2  shows  the  response  of  a 
deep  saturated  soil  column  to  two  field  earthquakes  (Loma  Prieta,  Hollister  Airport 
and  El  Centro),  one  artificial  time  history  (Folsom  Dam  Record  B)  and  a  centrifuge 
earthquake  base  input  motion  (Model  3c,  earthquake  1).  The  SHAKE  output  is 
plotted  in  terms  of  stress  ratio  against  depth,  and  compared  with 

a)  a  simple  linear  function  for  rd  (Iwasaki  et  al.  (1978), 
rd=l- 0.015  z  (z  in  metres), 

recommended  as  suitable  for  depths  up  to  25m,  and 

b)  a  log  function  proposed  based  on  this  study, 
rd  =  1 .25  -  In  z  (z  in  metres), 

which  fits  the  data  well  for  all  depths  greater  than  around  2m. 
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Figure  4.3  Shear  stress  ratio  required  at  different  depths  to  generate  10  -  20%  ground 
surface  acceleration,  based  on  rd  =  1  -  0.0 15z,  rd  =  1.25  -  0.2  In  (z)  and  Youd  and 
Idriss 

Figure  4.3  shows  a  comparison  between  these  two  functions  and  the  output  based  on  a 
more  complex  function  proposed  by  Youd  and  Idriss  (1997),  equation  (3),  for  ground' 
surface  accelerations  of  10%  and  20%g. 

The  true  character  of  stress  ratio  with  depth  is  significantly  affected  by  the 
development  of  excess  pore  pressure  in  the  soil  column.  In  practice,  the  ground 
surface  acceleration  time  history  will  be  distorted  by  the  liquefaction  in  the  soil  strata 
beneath,  which  often  isolates  the  surface  from  the  underlying  motion.  In  this  study 
the  base  input  motion  was  used,  and  scaled  using  an  amplification  factor  derived  from 
the  actual  profile  of  amplification  observed  in  the  specimen.  The  stress  history  at  any 
depth  was  then  computed  following  equation  (1). 

Following  this  approach,  the  full  dataset  was  processed.  Appendix  A  presents  the  data 
in  detail. 

4.4  SHAKE  analysis  of  soil  column 

SHAKE  may  also  be  used  for  the  prediction  of  time  histories  of  output  acceleration  at 
different  levels  in  the  soil  column.  However,  the  onset  of  liquefaction  has  a  dramatic 
effect  on  the  acceleration  in  the  upper  part  of  the  soil  column,  and  this  is  not  reflected 
in  the  SHAKE  calculation.  Fig.  4.4  shows  the  output  motion  for  layer  1  (surface) 
compared  with  the  actual  surface  motion  and  input  base  motion  for  Model  3c, 
earthquake  1  (see  also  Fig.  2.9). 
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SHAKE  output  (layer  1:  surface) 


Figure  4.4  Comparison  between  actual  surface  motion  and  SHAKE  calculation  in 
liquefied  specimen.  Model  3c,  eql 

The  SHAKE  soil  model  used  in  the  calculation  was  in  20  layers,  with  upper  bound 
values  of  soil  stiffness  and  damping.  The  soil  column  was  around  26  m  deep,  with  the 
water  table  at  the  ground  surface.  The  output  confirms  that  estimates  of  amplification 
close  to  1  are  realistic  for  the  purposes  of  the  analysis. 

4. 5  Comparison  with  published  Magnitude  Scaling  Factors 

Many  excess  pore  pressure  records  showed  a  cyclic  response  during  pore  pressure 
buildup,  and  the  peak  pore  pressure  would  reach  the  ultimate  limiting  value 
significantly  before  the  residual.  Typical  data  from  Model  3c,  earthquake  1  is  shown 
in  Figure  4.5.  The  time  at  which  the  residual  pore  pressure  reached  the  limiting  value 
was  adopted  as  the  most  consistent  and  reliable  definition  of  the  onset  of  initial 
liquefaction. 

Using  the  period  associated  with  the  dominant  frequency  of  shaking,  the  number  of 
‘cycles’  to  reach  liquefaction  was  calculated  based  on  the  time  and  the  dominant 
period  of  earth  shaking.  Table  4.1  illustrates  the  output. 
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Actual  depth 

Prototype 

Static  pore 

Depth  fluid 

Total 

Effective 

(mm) 

depth  (m) 

pressure  at 

on  surface 

vertical 

vertical 

50g  (KPa) 

(m) 

stress, 

stress,  a’v 

cv  (KPa) 

(KPa) 

439 

!  21.95 

263 

0.391 

457 

194 

Maximum 

Maximum 

Cyclic  stress 

Effective 

Onset  of 

Number  of 

cyclic  stress 

cyclic  stress 

ratio  (CSR) 

vertical 

initial 

cycles 

ratio  (raw) 

ratio 

x/a’v 

stress  (tsf) 

liquefaction 

equivalent 

(smoothed) 

(seconds) 

(N) 

0.232 

0.225 

0.146 

2.1 

0.308 

4.27 

Data  for : 

Model  3c 

PPT30 

G  level 

Unit  weight 

Unit  weight 

Fluid  unit 

Earthquake 

Frequency 

(dense) 

(loose) 

weight 

start 

(Hz) 

(seconds) 

50  g 

i  20.7  KN/m3 

20.2  KN/m3 

11.77 

0.15 

27 

KN/m3 

Table  4. 1  Base  data  and  calculations  of  stress  ratio  and  number  of  cycles  to 
liquefaction  for  PPT30,  Model  3c,  earthquake  1  (refer  also  to  Figure  4.4) 

The  base  shaking  motion  shown  in  Figure  4.5  was  selected  to  simplify  the  derivation 
of  a  representative  Ka  value.  The  selection  of  Magnitude  Scaling  Factors  was  based 
on  data  reproduced  by  Youd  and  Idriss  (1997).  It  was  found  that  the  original  Seed  and 
Idriss  MSF  values  generate  a  ‘flat’  curve  represented  by  the  function  F  =  5.24  N‘° 61. 
Proposals  by  Ambraseys  can  be  represented  by  a  relatively  steep  curve  matched  by 
the  function  F  =  3.48  e  “0  08N.  A  more  recent  proposal  by  Idriss  (Youd  and  Idriss, 
Table  3)  yields  a  steeper  curve,  represented  by  the  function  F  =  3.47  N  ~0,45. 
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Figure  4.5  Development  of  excess  pore  pressures  in  Model  3c,  earthquake  1 


Figures  4.6,  4.7  and  4.8  show  the  families  of  curves  generated  by  these  three  functions 
compared  with  the  centrifuge  data  of  normally  consolidated  50%  RD,  first  earthquake 
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specimens.  The  Idriss  MSF  factors  were  selected  as  most  appropriate  (of  the  three 
choices),  although  a  slightly  flatter  curve  could  have  been  chosen  to  match  the  data 
more  closely.  Using  the  Idriss  curve  function,  the  least  squares  best  fit  curve  to  the 
data  at  1  tsf  was  found,  and  this  is  shown  in  Figure  4.8. 


Shear  stress  ratio 


Figure  4.6  Comparison  between  all  data  and  the  shape  of  the  liquefaction  curves 
predicted  using  the  MSFs  of  Seed  &  Idriss  (data  only  reaching  90%  excess  pore 
pressure  or  greater) 

Shear  stress  ratio 


Number  of  equivalent  uniform  cycles  to  liquefaction 


Figure  4.7  Comparison  between  all  data  and  the  shape  of  the  liquefaction  curves 
predicted  using  the  MSFs  by  Ambraseys 
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Shear  stress  ratio 


Number  of  equivalent  uniform  cycles  to  liquefaction 

Figure  4.8  Comparison  between  first  earthquake  data  and  the  shape  of  the 
liquefaction  curves  predicted  using  the  MSFs  by  Idriss 

Using  the  same  family  of  curves,  each  data  point  was  then  scaled  to  its  equivalent 
stress  ratio  at  10  cycles.  As 

(t/ct’v)n  =  (t/o’v)n=15  .  A  N'b 

where  A  and  B  are  the  factors  describing  the  family  of  curves,  and  therefore 

(t/ct’v)n=10  =  (t/ct’v)n=15  •  A  10'B 

then 

(t/ct’v)n=io  =  (t/c’v)n  .  10'b/N'b 


4. 6  Comparison  between  centrifuge  Kg  and  laboratory  Kg  results 

The  ratio  of  the  equivalent  cyclic  stress  ratio  at  ten  cycles  (defined  as  the  liquefaction 
cyclic  resistance  strength  ratio  CRR)  to  the  cyclic  stress  ratio  at  10  cycles  for  the  best 
fit  ltsf  curve  was  calculated,  giving  K^,  as  shown  in  Figure  4.9. 
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Kct 


Vertical  effective  stress,  tsf 

Figure  4.9  Comparison  between  first  earthquake  data  (50%  RD,  normally 
consolidated  only),  expressed  as  Ko  at  10  cycles,  compared  with  Olsen  (1984)  and 
laboratory  test  data  from  Vaid  for  Tailings  and  Ottawa  sand 

Figure  4.9  also  includes  data  of  calculated  by  Hynes  and  Olsen  from  laboratory 
test  data  of  Vaid  on  Tailings  sand  and  Ottawa  sand.  The  curves  shown  in  the  figure 
were  calculated  using  Olsen’s  function  KCT  =  (cr’y/'1  ,  Olsen  (1984). 

Shear  stress  ratio 


Number  of  equivalent  uniform  cycles  to  reach  50%  excess  pore  pressure 

Figure  4.10  Comparison  between  all  data  and  the  shape  of  the  liquefaction  curves 
predicted  using  the  MSFs  of  Seed  &  Idriss  at  50%  excess  pore  pressure,  first 
earthquakes  only 
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Shear  stress  ratio 


Number  of  equivalent  uniform  cycles  to  reach  50%  excess  pore  pressure 

Figure  4. 1 1  Comparison  between  all  first  earthquake  data  and  the  shape  of  the 
liquefaction  curves  predicted  using  the  MSFs  by  Idriss  at  50%  excess  pore  pressure 

Ko 


Figure  4.12  Comparison  between  first  earthquake  data  (50%  RD,  normally 
consolidated  only)  expressed  as  Kcr  at  10  cycles  for  50%  excess  pore  pressure, 
compared  with  Olsen  (1984) ,  calculated  using  Seed  MSFs 
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The  range  of  data  which  showed  90%  plus  excess  pore  pressure  generation  was 
limited  to  relatively  shallow  depths.  A  second  calculation  was  made  to  assess  the 
value  of  Ka  at  50%  excess  pore  pressure  development.  In  this  case,  although  the  data 
are  more  constrained  in  terms  of  number  of  cycles  to  reach  50%,  it  is  clear  from 
Figures  4.10  and  1 1  that  neither  the  Idriss  MSF  or  the  Seed  and  Idriss  equivalent 
follow  the  trend  of  the  data  at  both  low  and  high  numbers  of  cycles.  This  has  a 
significant  effect  on  the  calculation  of  Ka,  as  shown  in  Figures  4.12  and  4.13. 

Ka 


Figure  4.13  Comparison  between  first  earthquake  data  (50%  RD,  normally 
consolidated  only)  expressed  as  Ka  at  10  cycles  for  50%  excess  pore  pressure, 
compared  with  Olsen  (1984) ,  calculated  using  Idriss  MSFs 

It  is  clear  that  the  values  of  Ka  are  highly  sensitive  to  the  selection  of  the  shape  of  the 
characteristic  curve,  which  is  markedly  flatter  at  high  numbers  of  cycles  than  the 
standard  MSFs  would  predict.  The  variability  of  Ka  is  widely  remarked  upon  in  the 
literature,  Youd  and  Idriss  (1997). 

4. 7  Comparison  with  ‘Stress  focus  ’  concept 

Alternative  representations  of  the  CRR  and  CR  as  a  function  of  depth  have  been 
proposed  in  order  to  support  the  extrapolation  of  Ks  from  low  to  high  effective 
confining  pressures.  Hynes  (1998)  describe  the  use  of  log-log  plots  of  CRR  and  CR 
vs  depth  to  indicate  the  'stress  focus',  towards  which  liquefaction  resistance  measured 
under  low  effective  confining  stress  may  be  extrapolated  to  deduce  resistance  under 
high  effective  confining  stress. 
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Figure  4. 14  Stress  focus  plot  using  centrifuge  data,  normally  consolidated,  50%  RD, 
computed  for  1 0  cycles  using  Idriss  MSF,  and  compared  with  VELACS  laboratory 
test  data  on  Nevada  sand  corrected  to  50%  RD  (ClUCyclic  squares,  CSS  diamonds) 
and  projected  to  10  cycles  using  Seed  MSF 

The  data  from  the  centrifuge  experiments  are  shown  in  log  CRR  -  log  depth  form  in 
Figure  4.14.  Two  curves,  representing  f  =  0.7  and  f  =  0.9  were  drawn  through  the 
CRR  at  1  tsf.  (For  the  centrifuge  data,  this  was  the  best  fit  value.  For  Vaid’s  data,  this 
was  the  1  tsf  value.)  Also  shown  in  the  figure  are  laboratory  test  data  carried  out  by 
Earth  Technology  Corporation  for  the  VELACS  project. 

[  The  laboratory  test  data  were  processed  as  follows.  Seven  experiments  were 
identified  from  the  VELACS  data,  three  undrained  cyclic  triaxial  tests  (ClUCyclic) 
and  four  undrained  cyclic  simple  shear  tests  (CSS).  These  were  conducted  on  sand  at 
60%  RD.  The  triaxial  data  was  converted  to  a  simple  shear  equivalent  using  a  factor 
of  2/root3  x  q/2  /a^,  following  Castro.  The  data  was  then  further  adjusted  for  relative 
density  using  the  Holtz  and  Gibbs  relation  for  SPT  and  relative  density  for  coarse 
sand.  At  1  tsf,  the  number  of  blows  at  50%  RD  is  around  10.  At  60%,  this  rises  to 
around  13.4.  From  Seed  for  a  Magnitude  7.5  earthquake  the  CSR  at  50%  is  around 
0.102  and  at  60%  is  around  0.14.  hence  the  data  is  scaled  by  a  factor  of  0.729. 

Finally,  the  data  was  extrapolated  to  the  equivalent  CSR  at  10  cycles.  This  was  done 
using  both  the  Idriss  curve  and  the  flatter  Seed  curve.  It  was  found  that  the  flatter 
Seed  curves  collapsed  the  laboratory  data  on  top  of  the  centrifuge  data  much  more 
closely  than  the  Idriss  curve  (which  more  closely  reflected  the  centrifuge  data).  These 
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laboratory  test  data  for  Nevada  sand  were  then  plotted  in  Figure  4.15,  using  solid 
squares  for  the  triaxial  data  and  open  diamonds  for  the  simple  shear  data.  ] 

Liquefaction  Cyclic  Resistance  Strength  Ratio  (CRR) 

0.01  0.1  1 


Figure  4.15  Stress  focus  plot  using  centrifuge  data,  normally  consolidated,  50%  RD, 
computed  using  Idriss  MSF  for  10  cycles,  and  compared  with  VELACS  laboratory 
test  data  on  Nevada  sand  corrected  to  50%  RD  (ClUCyclic  squares,  CSS  diamonds) 
and  projected  to  10  cycles  using  Seed  MSF 

Figure  4.15  shows  the  comparable  plot  of  CR  as  a  function  of  depth.  Recalling  that 
the  CRR  and  CR  are  calculated  essentially  by  dividing  the  applied  cyclic  shear  stress 
at  initial  liquefaction  by  the  initial  effective  confining  stress  (and  scaling  to  ten 
cycles),  it  is  instructive  to  compare  the  data  with  'true'  stress  ratios  denoted  as 
mobilised  angles  of  shearing  resistance  <{>'  in  the  conventional  manner.  These  may  be 
simply  added  as  straight  lines  in  both  figures. 


4.8  Stress  path  analysis 

A  full  derivation  of  the  stress  path  analysis  is  presented  in  Appendix  B.  The  stress 
path  for  a  typical  soil  element  is  shown  in  Figure  4.16,  for  Model  3c,  earthquake  1, 
depth  22m  (Table  3.2).  The  value  of  Ko  based  on  Ko  =  1  -  sin  <j>',  where  (j)'  =  33 
degrees  as  the  boundary  conditions  for  a  soil  element  in  the  centrifuge  specimen  are 
the  equivalent  of  zero  lateral  strain.  Separate  plots  show  the  mean  effective  confining 
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stress,  s',  the  shear  stress  t  and  applied  cyclic  shear  stress  xav,  and  the  'true'  mobilised 
angle  of  shearing  resistance.  The  angle  of  shearing  resistance  computed  for  this  data 
point  based  on  the  liquefaction  cyclic  resistance  ratio,  with  CRR  =  0.0868,  CR  = 
0.157,  avo'  -  194  KPa  (1.81  tsf ),  is  given  by  sin  4>’  =  CRR/(1-CRR),  or  <(>'  =  5.45 
degrees. 

This  data  point  reached  the  origin  and  a  near  zero  effective  stress  state.  However 
others,  at  higher  effective  confining  stress,  were  observed  not  to  reach  the  same  state 
and  despite  continued  cycling,  the  generation  of  excess  pore  pressure  was  limited. 

One  such  data-point  was  Model  4d,  earthquake  2,  shown  in  Figure  4.1 7.  Here  the 
stress  path  moves  progressively  towards  the  origin,  excess  pore  pressure  developing  at 
a  steady  rate.  After  a  number  of  cycles,  excess  pore  pressure  stops  rising  and  the 
downward  migration  of  the  stress  path  is  arrested.  This  point  reached  around  44%  of 
the  initial  vertical  effective  stress. 
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Figure  4.16  Stress  path  for  soil  element  at  initial  vertical  effective  stress  of  194  KPa 
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4. 9  Data  of  excess  pore  pressures 

Under  high  effective  confining  stresses,  it  was  noted  that  in  many  cases,  excess  pore 
pressures  did  not  reach  the  initial  vertical  effective  stress,  as  was  the  case  with  data 
from  shallow  depths.  Plotting  all  data  points  in  Figure  4.18,  it  is  clear  that  there  is  a 
reducing  trend  at  initial  vertical  effective  stresses  above  around  300  KPa  and  excess 
pore  pressures  do  not  reach  sufficient  levels  to  cause  initial  liquefaction,  despite 
continued  cycling.  In  this  figure,  data  points  include  data  of  all  densities,  first  and 
following  earthquakes,  normally  and  overconsolidated  specimens. 
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Figure  4.18  Comparison  between  upper  limit  of  excess  pore  pressure  generation  and 
initil  effective  vertical  stress,  all  data 

A  similar  phenomenon  has  been  observed  in  torsional  simple  shear  tests,  Ishihara 
( 1 996)  from  Ishihara  and  Nagase  (1988).  Fig.  1 9  shows  that  although  there  is  a 
minor  effect  of  overconsolidation,  which  appears  to  reduce  the  magnitude  of  the 
excess  pore  pressure  slightly  below  that  of  the  normally  consolidated  data,  as  might 
be  expected.  The  trend  with  increasing  effective  vertical  stress  is  identical,  however. 
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avo' 


Figure  4.19  Comparison  between  upper  limit  of  excess  pore  pressure  generation  and 
initial  effective  vertical  stress,  first  earthquake  data,  normally  and  overconsolidated 


tfvo' 

Figure  4.20  Ru  as  a  function  of  stress  ratio,  for  first  earthquake  <50%  RD,  normally 
consolidated  data 

Figure  4.20  shows  a  selection  of  the  data,  concentrating  only  on  normally 
consolidated  data  from  first  earthquakes  at  50%  Relative  Density  or  less.  At  an  initial 
vertical  effective  stress  of  less  than  around  300  KPa,  even  low  applied  cyclic  stress 
ratios  are  sufficient  to  achieve  excess  pore  pressures  of  100%.  At  higher  vertical 
effective  stresses,  however,  this  is  not  the  case  and  Ru  falls  progressively.  Although 
most  of  the  data  is  around  50%  RD,  the  very  highest  effective  confining  stresses  are 
associated  with  a  low  relative  density  (39%),  which  may  explain  the  higher  values  of 
Ru  than  suggested  by  the  general  trend. 
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This  phenomenon  is  not  predicted  using  standard  design  methods. 

Fig.  4.21  breaks  the  data  out  by  earthquake,  from  which  it  is  clear  that  there  is  little  to 
distinguish  between  first  and  succeeding  shaking  events.  The  implication  of  this  is 
that  the  pattern  of  excess  pore  pressure  generation  is  very  similar  in  successive 
earthquakes,  showing  the  same  limit  and  trend  with  depth. 


Figure  4.21  Comparison  between  upper  limit  of  excess  pore  pressure  generation  and 
initial  effective  vertical  stress,  overconsolidated,  by  earthquake 


Figure  4.22  CSR  as  function  of  cycles  and  Ru  (%),  first  earthquakes,  <50%  RD, 
normally  consolidated  data,  100  KPa  (~1  tsf) 
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Finally,  the  relationship  between  stres  ratio  and  number  of  cycles  to  reach  different 
levels  of  excess  pore  pressure,  expressed  as  Ru  (%)  is  shown  in  Fig.  4.22  .  The  curves 
appear  to  be  broadly  parallel  to  each  other,  rather  than  tending  to  merge  at  low  values 
of  CSR  as  would  be  predicted  from  laboratory  element  tests.  This  confirms  that  the 
observation  of  limited  Ru  with  depth  is  not  simply  a  consequence  of  low  CSR. 
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5  CONCLUSIONS 

1 .  A  large  database  of  the  development  of  excess  pore  pressure  in  a  uniform 
saturated  sand  bed  has  been  established. 

2.  The  extrapolation  of  the  simplified  procedure  for  the  evaluation  of  liquefaction 
resistance  using  the  correction  factor  Ko  may  not  be  warranted  at  high  confining 
stresses.  The  centrifuge  derived  values  of  Ka  do  not  match  the  laboratory 
developed  values.  This  is  probably  because  the  laboratory  cyclic  triaxial  test  does 
not  provide  similar  boundary  conditions  to  the  centrifuge  model,  which  are 
considered  to  be  more  representative  of  the  field  conditions. 

3.  Observation  of  the  stress  path  followed  by  soil  elements  as  they  approach  a 
condition  of  near  zero  effective  stress  shows  that  the  strength  of  the  soil  is  not 
related  to  the  initial  effective  vertical  stress. 

4.  There  is  a  cut-off  in  the  development  of  excess  pore  pressure  which  is  a  function 
of  the  initial  vertical  effective  stress.  Further  studies  are  required  to  confirm  that 
this  cap  is  likely  also  to  be  a  function  of  earthquake  amplitude  and  is  probably 
associated  with  strain  amplitude. 

5.  It  is  postulated  that  the  relationship  between  Ru  and  shear  strain  derived  from 
laboratory  element  test  data  may  be  conservative  when  applied  to  centrifuge  or 
field  boundary  conditions. 

6.  The  implications  of  this  finding  are  significant  for  the  assessment  of  the  risk  posed 
by  earthquakes  to  large  earth  dams  and  similar  structures,  where  the  high  effective 
confining  stresses  in  the  foundations  may  lead  to  an  overly  conservative 
assessment  of  the  requirement  for  remediation. 
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Stress  path  analysis 
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APPENDIX  B:  CALCULATION  OF  THE  STRESS  PATH 

The  stress  path  followed  by  an  element  in  the  sand  bed  was  calculated  in  s’,  t  space 
using  the  time  history  of  shear  stress  computed  from  the  acceleration  record.  Two 
models  were  selected  for  first  analysis,  Model  3c,  earthquake  1  and  Model  5d, 
earthquake  1 .  In  Model  3c,  ppt  30  showed  clear  evidence  of  liquefaction  after  4.3 
cycles,  reaching  100%  of  the  initial  vertical  effective  stress  with  a  cyclic  stress  ratio 
t/ovo’  =  0.155.  As  described  previously,  Steedman  (1999a),  the  cyclic  stress  ratio  is 
calculated  as  0.65  times  the  peak  stress  ratio  of  0.239,  following  conventional 
practice.  In  the  case  of  Model  3c,  the  peak  stress  ratio  was  based  on  the  smoothed 
acceleration  input  from  acc  12609  (at  the  bottom  of  the  sand  bed).  In  Model  5d,  the 
pore  pressure  in  the  loose  layer  reached  a  peak  level  of  around  52%  of  the  initial 
effective  vertical  stress,  partway  through  the  shaking,  and  then  was  capped.  This 
observation  prompted  a  closer  examination  of  the  stress  paths  being  followed  at 
these  locations  in  the  models. 


Figure  B.1 


To  calculate  the  stress  path,  the  stress  ratio  of  shear  stress  to  initial  effective  vertical 
stress  was  used  first  to  compute  the  applied  shear  stress  on  a  horizontal  plane  as  a 
function  of  time,  using  the  time  history.  The  stress  on  a  vertical  plane  is  given  by 

oh’  =  Ko  av’. 

Shear  stress  applied  to  horizontal  (and  vertical)  planes  causes  the  Mohr  circle  to 
dilate,  without  change  in  the  mean  effective  stress,  s’,  as  shown  in  Figure  B.1. 

In  this  (elastic)  condition,  prior  to  yield,  it  may  be  shown  that 

s’  =  0.5  ov’  (1  +  Ko),  and 

t  =  0.5 { av’ (1  -Ko)2  +4t2}05 

Once  the  stress  state  reaches  a  condition  of  yield,  denoted  by  a  maximum  mobilised 
angle  of  shearing  resistance,  <j)'max,  as  shown  in  Figure  B.2,  then  further  increase  in 
the  applied  shear  stress  on  a  horizontal  plane  forces  the  mean  effective  confining 


B-l 
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stress  to  increase,  as  the  circle  is  constrained  by  the  <J>’max  line.  In  this  ‘yielding’ 
condition  it  may  be  shown  that 

s’  =  crv’{sec2  <t>’max  ±  sec  <j>’max  V(tan2  <{>’max  -  tan2  p )} ,  and 

t  =  s’  sin  fmax 

where  tan  p  =  x/av’  =  (t/gvo’)  avo7av’ .  A  more  extensive  discussion  of  this  approach 
as  adopted  for  the  calculation  of  stress  paths  behind  a  retaining  wall  was  presented 
by  Steedman  (1999b). 
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Figure  B.2  Example  of  the  dilation  of  a  ‘yielding’  stress  state 

At  any  stage,  the  mobilised  angle  of  shearing  resistance  may  be  simply  calculated  as 

<{)’mob  =  asin(s’/t) 

As  excess  pore  pressure  develops,  the  vertical  effective  stress  is  reduced,  and  s’  is 
also  reduced.  Assuming  that  Kq  is  constant,  then  the  radius  of  the  Mohr’s  circle, 
given  by  t,  must  also  reduce.  Depending  on  the  amplitude  of  the  applied  shear 
stress  and  the  rate  of  generation  of  excess  pore  pressure,  the  stress  path  may  be 
driven  onto  the  yield  surface  or  not. 

The  parameter  Ko  is  seen  to  be  fundamental  to  the  prediction  of  the  stress  path,  both 
in  its  initial  value  and  in  the  assumption  that  it  is  constant  at  any  stage  prior  to 
liquefaction. 

Following  this  approach,  the  stress  paths  were  computed  for  Model  3c  and  5d. 

Figure  B.3  shows  the  pore  pressure,  acceleration  and  stress  ratio  for  Model  3c 
associated  with  ppt  30,  at  mid  depth  in  the  loose  layer,  under  an  initial  vertical 
effective  stress  of  194  KPa  (1.8  tsf). 

Two  values  of  Kq  were  used,  Ko  =  1  and  K0  =  1  -  sin  <j>’max .  Figure  B.4  shows  the 
stress  path  in  s't  space  for  ppt  30,  assuming  Ko  =  1.  Figures  B.5  and  B.6  show  the 
progression  over  time  of  s’  and  t  respectively.  Figure  B.7  shows  the  mobilised  angle 
of  shearing  resistance. 
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Figure  B.3  Time  histories  of  pore  pressure,  acceleration  and  shear  stress  ratio, 
Model  3c 


Figure  B.4  Stress  path  at  mid-depth  in  loose  layer,  Model  3c,  earthquake  1  (Ko  =1) 
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Figure  B.7  Mobilised  angle  of  shearing  resistance,  Model  3c  (Ko=1) 

For  Ko  =  1-sin  <j)’max,  then  Ko  =  0.455.  This  has  a  dramatic  effect  on  the  predicted 
stress  path  in  s’,  t  space,  Figure  B.8. 


Figure  B.8  Stress  path  for  Ko  =  0.455,  Model  3c,  earthquake  1 


Figure  B.9  Mean  effective  confining  stress,  s’,  for  Ko  =  0.455,  Model 


Figure  B.10  Shear  stress  for  Ko  =  0.455,  Model  3c 
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Figure  B.1 1  Mobilised  angle  of  shearing  resistance,  Ko  =  0.455,  Model  3c 


Model  5d  showed  a  rise  in  excess  pore  pressure  in  the  loose  layer,  which  reached  a 
limiting  value  partway  through  the  shaking  event.  Such  behaviour  has  been  oberved 
in  laboratory  tests,  Ishihara  (1996). 

Figure  B.1 2  shows  the  pore  pressure  at  mid-depth  in  the  loose  layer,  the 
acceleration  and  shear  stress  time  histories  for  Model  5d,  earthquake  1.  The  initial 
vertical  effective  stress  was  1037  KPa  (9.7  tsf). 
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Figure  B.12  Time  histories  of  pore  pressure,  acceleration  and  applied  shear  stress 
for  Model  5d,  earthquake  1 ,  at  mid-depth  in  the  loose  layer 

Figure  B.12  also  shows  a  comparison  between  the  raw  and  filtered  (Butterworth 
band  pass  filter,  low  250  Hz,  high  10  Hz)  acceleration  records.  The  filtered  time 
history  was  used  for  computing  the  applied  shear  stress  because  it  was  considered 
to  be  a  more  accurate  representation  of  the  true  strain  history  in  the  soil  specimen, 
without  the  high  frequency  noise  distortion,  particularly  towards  the  end  of  the 
earthquake. 
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For  Model  5d,  the  analysis  shows  that  with  the  stress  path  does  not  reach  a  yield 
condition,  staying  well  below  the  yield  surface  for  Ko  =  1,  Figure  B.13. 
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Figure  B.13  Stress  path  for  Model  5d,  earthquake  1,  Ko  =  1 
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Figure  B.14  mean  effective  confining  stress,  Model  5d,  earthquake  1,  Ko  =  1 


Ko  =1,  <j>’  =  33, 

ovD  =  1280  KPa 
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Figure  B.15  Shear  stress  t  and  applied  shear  stress  t  for  Model  5d,  earthquake  1 , 
Ko  =  1 


20 


Figure  B.16  Mobilised  angle  of  shearing  resistance,  Model  5d,  earthquake  1,  Ko  =  1 

The  detail  of  the  stress  path  was  examined  by  windowing  the  data  to  examine  a 
short  period  during  the  shaking  when  the  mean  effective  stress  was  ‘capped’  at 
around  52%.  The  stress  path  for  around  two  cycles  of  s’  is  shown  in  Figure  B.17. 
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1150 

Figure  B.17  Stress  path  for  two  cycles  of  s’,  Model  5d,  earthquake  1,  Ko  =  1 

Also  shown  in  Figure  B.17  is  a  theoretical  shape  for  the  path,  assuming  that  s’  and  t 
may  be  represented  by  sine  functions.  Figure  B.18  and  B.19  show  the  paths  in 
mean  effective  stress  and  shear  stress,  together  with  the  arbitrary  sine  functions 
used  to  model  the  fluctuations  with  time.  In  all  of  these  figures,  the  x  axis  is  the 
number  of  data  points. 
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Figure  B.18  Enlargement  showing  detail  of  s’  and  t  during  ‘capped’  portion  of 
shaking 
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Figure  B.19  Comparison  between  applied  shear  stress  and  t 


Figure  B.20  Mobilised  angle  of  shearing  resistance 

Figure  B.20  shows  the  mobilised  angle  of  shearing  resistance  for  Model  5d, 
earthquake  1,  for  the  condition  Ko  =  1. 
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Similarly  to  Model  3c,  the  stress  path  at  mid-depth  in  the  loose  layer  is  changed 
dramatically  in  the  case  of  a  smaller  Ko.  Using  a  new  value  of  Ko  =  0.455,  the  stress 
path  was  recalculated  for  Model  5d,  earthquake  1,  Figure  B.21. 


Figure  B.21  Stress  path  for  Model  5d,  earthquake  1 ,  for  Ko  =  0.455 


Figure  B.22  shows  the  shear  stress  as  a  function  of  time. 


B-12 


Report  RSS/J9901R08  Rev  0 


400 
300 
200 
t,x  100 
0 

-100 
-200 

1  201  401  601  801  1001  1201  1401  1601  1801  2001 

Figure  B.23  Comparison  between  applied  shear  stress  and  t  for  Model  5d, 
earthquake  1,  Ko  =  0.455 
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Figure  B.24  Mobilised  angle  of  shearing  resistance,  Model  5d,  earthquake  1 ,  Ko  = 
0.455 
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